Ageing of a granular pile induced by thermal cycling 
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Abstract. Here we show that variations of temperature, even of a few degrees in amplitude, induce the ageing of a granular 
pile. In particular, we report measurements of physical properties of a granular heap submitted to thermal cycles. Namely, we 
focus on the evolution of the thermal linear-expansion coefficient and of the thermal conductivity of the pile with the number 
of cycles. The present contribution nicely supplements a recent article we published elsewhere [Phys. Rev. Lett. 101, 148303 
(2008)] and introduces a different and promising method to impose temperature cycles to a granular pile. 
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INTRODUCTION 

A static granular pile is essentially a fragile construc- 
tion Oil, |2fl; Minute perturbations applied at the micro- 
scopic scale (the surface roughness of a grain) can lead 
to macroscopic reorganisations of the pile [3]. A remark- 
able signature of this fragility is the sensitivity of a gran- 
ular pile to temperature variations. First presented as a 
hindrance to assess the sound propagation in sand J3l or 
to perform reproducible measurements of the stress dis- 
tribution at the base of a granular pile JH, the sensitivity 
of granular matter to temperature fluctuations has been 
since used to pack grains [7, 8, 9]. Indeed, a granular 
heap submitted to temperature cycles, even of a few de- 
grees in amplitude HQ ESDI, experiences successive 
large-scale "static avalanches" which induces the 
slow compaction of this fragile construction. Such reor- 
ganisations are possible because of the stress anisotropy 
inside the pile IU2I1 and of the surface properties of the 
grains 

In a recent letter [11], we reported a time-resolved 
study of the dynamics associated with the compaction of 
a granular column submitted to thermal cycles. Here, we 
demonstrate that temperature variations, even of a few 
degrees in amplitude, also induce ageing in a granular 
pile exactly as moisture 111311 . constant applied-stress lfl4ll 
and chemical reactions between grains IU5I1 do. Here, we 
dwell on two physical properties of the granular pile, 
namely Kg, its thermal linear-expansion coefficient, and 
X g , its thermal conductivity which both evolve in time 
because of temperature varations. Both physical prop- 
erties remain poorly studied in the case of a 3D pile 
free to dilate and, thus, to reorganise under tempera- 
ture variations. Indeed, previously, only X g has been as- 
sessed numerically and experimentaly in the case of a 
bidimensional and compressed static-bed (see 11611 and 
references therein). We thus address here the issue of the 



dependence of JQ, and Xb with the numbers of imposed 
temperature-cycles. 
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FIGURE 1. Sketch of the first experimental setup. Inset: 
picture of the upper part of the column. The granular level is 
indicated by the white dotted-line. 



THERMAL DILATION 

In this section, we report, in a first experimental configu- 
ration, the evolution of the thermal linear-expansion co- 
efficient of a granular column submitted to thermal cy- 
cles. The experiment consists in imposing temperature 
cycles to a granular column and in measuring the result- 
ing variations of the column height. 



Experimental setup 

The experimental set-up (Fig. Q]) consists in a ver- 
tical glass tube (height 1.7 m, inner diameter 13 mm) 
firmly fastened to a wall into the basement of the physics 
department, in order to flee mechanical vibrations. The 
sample consists of a column (height H) of spherical glass 
beads (diameter d) poured into the tube. An additional 
gas-input, at the bottom of the column, makes possible 
to loosen the pile thanks to an upward flow of dry ni- 
trogen. The temperature cycles are imposed by means 
of a heating cable (Prolabo, 40 W/m) directly taped on 
the outer surface of the tube wall. The resulting temper- 
ature is measured by means of a sensor (PtlOO, located 
close to the free surface of the granular material) and a 
multimeter (Keithley, 196). The free surface of the ma- 
terial, which is illuminated by a red LED (Kingbright, 
L-793SRC-E, located inside the tube, above the granu- 
lar material) is imaged from the side with a video cam- 
era (Panasonic, WV-BP500) connected to a frame grab- 
ber board (Data Translation, DT2255). A macro, running 
under a data-processing software (WaveMetrics, IGOR 
Pro 4.0), drives the heating power, records the result- 
ing variations of the temperature and measures accu- 
rately the height H from the images: A subpixel resolu- 
tion (namely, less than a tenth of a pixel which typically 
stands for 5 jJ.m) is achieved by considering the average 
position of the free surface, assumed to correspond to the 
inflection point in the vertical intensity-profile averaged 
over the whole diameter of the tube. Measurements are 
performed 20 times per temperature cycle. 

Due to long experimental times, we limit our report to 
a given diameter d = (5 1 ± 90) /xm of the grains (Matra- 
sur Corp.) and to a given period 2n/co = 600 s of the cy- 
cles. The cycling period, 10 minutes, is arbitrarily chosen 
to be small enough to avoid excessively-long experimen- 
tal times but large enough to insure that the associated 
thermal penetration-length l p = \J2X /(Ceo) ~ 6 mm is 
about the tube radius (A ~ 0.2 W m" 1 K" 1 and C ~ 10 6 J 
m~ 3 K _1 respectively denote the thermal conductivity 
and heat capacity of a typical glass-grains pile lfl7ll .) It 
is here crucial to note that the column is heated homoge- 
neously along its whole length but that the temperature 
is likely to vary in the radial direction. 

Prior to each experiment, the granular column is pre- 
pared in a low-density state thanks to the dry-nitrogen 
upward flow. The top of the column is then higher than 
the field imaged by the camera (typically 1 cm above) 
and we set the amplitude of the cycles, AT, to the largest 
accessible value, AT = 27.1° C. The preparation of the 
sample ends when the top of the column enters the ob- 
servation field. At this point, the granular column is 
"quenched": The amplitude of the cycles is set to the 
chosen value Ar lying between 10.8 and 27.1°C, which 
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FIGURE 2. Height variation h„ vs. number of cycles 

n. One observes first an exponential behavior at short time 
followed by a subsequent logarithmic creep at long time 
[The black curve corresponds to the test function h' n = ho + 
h e exp(— n/tic) + hi ln(«).] Inset - Oscillations of the column 
height associated with the temperature cycles : A„ and 8„ are 
respectively defined to be the amplitude of the increase and 
the drift of h„ at the cycle n (H = 140 cm, 2n/w = 600 s and 
AT= 10.8°C.) 



defines the origin of time t = 0. The granular column is 
subsequently submitted to, at least, 1000 cycles (7 days). 



Results 

Under the action of the temperature cycles the 
column height decreases: we report the variation 
h„ = H(2n n/co) — H(Q), where n denotes the time in 
units of the cycle period or, equivalently when integer, 
the number of imposed cycles (Fig. We observe 
that the thermal-induced compaction is a very slow 
phenomenon: after 7 days (1000 cycles), the decrease of 
the height is of about 1 .5 cm (about 1 % of the height H), 
which indicates that the system remains very far from 
the maximum compaction (roughly a decrease of about 
10% of the column height, i.e. h„ ~ 10 cm) within the 
experimental time. Accordingly, after the experiment, 
we checked that a single finger tap produces a collapse 
of the granular height of one centimeter at a rough esti- 
mate. We also checked during 3 days that the height H 
of the column stays constant when no temperature cycles 
are imposed, which proves that ambient mechanical 
vibrations and changes in the room temperature have no 
(or little) effect in our experimental conditions. 

Interestingly, the measurements are accurate enough 
to reveal the oscillations of H associated with the tem- 
perature variations (Fig. [2] inset). We observe on the raw 
data that the amplitude A n , which is proportional to Ar, 
increases logarithmically with n (Fig. [3]). The oscillations 
of H are due to the thermal dilation of both the tube and 




the granular material. In order to assess the contribution 
of the granular material alone, we first determine the am- 
plitude, 8h t (z), of the tube displacement in the laboratory 
frame as a function of the height z (origin at the bottom 
of the column) by marking its outer wall. The amplitude 
8h t (z) is found to be linear in z and the slope provides us 
with an estimate of the linear thermal-expansion coeffi- 
cient, K = (3.6 ±0.4) x 1CT 6 K -1 , of the tube material. 
Then, considering the relative variation of the inner vol- 
ume, we write the relation between the amplitude A„ and 
the relative variation, 8V g /V g , of the volume V g of the 
granular material: A„ - 8h,{0) = H(8V g /V g - 2kAT). 
Experiments performed for different height H demon- 
strate that A„ — 8h t (0) is proportional to H, which shows 
that 8Vg/V g is independent of H and, thus, that the whole 
height H of the granular column is involved in the ob- 
served oscillations of the free surface. For a homoge- 
neous temperature-variation, AT, in the whole cross- 
section of the tube, we would have 8V g /V g = 3K g AT, 
With this assumption, we obtain a rough estimate of the 
thermal expansion coefficient of the granular material, 
Kg ~ [3.4 + 0.03 ln(n)] x 10~ 6 K _1 , which is thus found 
to be very close to K and to increase logarithmically by 
about 5% during the first 1000 cycles (Pig. [3j as a conse- 
quence of a variation of about 1% of the density (Pig. [2j. 

THERMAL CONDUCTIVITY 

In this section, using a second experimental setup, we 
estimate the thermal conductivity of a granular pile and 
perform a time resolved study of its evolution under cy- 
cles of temperature. The experiment consists in heating 
and in measuring the temperature of a granular column 
along its axis of revolution thanks to a thin heating wire. 



FIGURE 4. Sketch of the second experimental setup. 

Experimental setup 

The experimental setup (Fig. HJ consists in a cop- 
per tube (inner diameter 1 cm; outer diameter 5 cm; 
length 15 cm) partially filled with glass beads (diameter 
d G [250 — 425] fim). The tube is maintained at a con- 
stant temperature T e = 55°C (numerical PID controller 
written in C++, precision 0.1 °C) whereas the granular 
column is heated along its axis thanks to a nickel wire 
(diameter r w =100 fim) connected to a power supply 
(Hewlett Packard, 6633A). We measure simultaneously 
the imposed current / and the resulting voltage V (Keith- 
ley, 196) in order to deduce the heating power P = UI 
and the wire temperature T w which we deduce from the 
resistance R w = U /I. The thermal conductivity X g of the 
granular material is then assessed from the temperature 
difference T w — T e which depends linearly on P/X g B17I1 . 



Results 

In order to measure the thermal conductivity, P is 
increased by steps (typically 10) up to T w — T e ~ 40°C 
and then decreased by steps down to 0. For each value of 
P, we wait until the steady state is reached and measure 
the associated voltage U and current /. We obtain X g = 
(0.16 ±0.02) W/m/K. This value is compatible with the 
value expected for a pile of glass beads {X g i ass =1.4 
W/m/K) surrounded by air (X air = 0.025 W/m/K) d. 

We interestingly also observed that, for the same ma- 
terial, X g slightly depends on the preparation (Fig. |5]i : 
X g ~ 0. 162 W/m/K if the system is tapped previous to the 
measurement, X g ~ 0.156 W/m/K if not. It is then partic- 
ulary interesting to consider the behaviour of the sample 
when subjected to several temperature cycles. When the 
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FIGURE 5. Thermal conductivity X g vs. number of cycles 

n. The upper curve corresponds to a pile which has been gen- 
tly tapped previous to the experiment. The lower curve corre- 
sponds to a loose pile. Note first that the conductivity is larger 
for a larger density. Moreover, the conductivity of the loose 
material increases significantly when one imposes the thermal 
cycles ( d e [250-425], Ar = 40°C.) 

measurements are repeated several times, one observes 
that the thermal conductivity of the loose sample signif- 
icantly increases with the number n of imposed cycles. 
By contrast, the conductivity of the tapped sample only 
slightly fluctuates around a constant value. 

The second experimental configuration does not make 
possible to observe the grains. However, a few experi- 
ments reported in details in I7J proved that the granular 
column is subjected to compaction even if the container 
does not dilate (heating cable along the column axis). 
Our results show that the increase in the thermal con- 
ductivity kg with the number of cycles n originates in the 
slow compaction of the material. 

CONCLUSION AND OUTLOOKS 

We reported measurements of the thermal linear- 
expansion coefficient Kg and of the thermal conductivity 
X g of a granular pile. We have shown that temperature 
variations, even of a few degrees in amplitude, can lead 
to the compaction of the pile and, as consequences, to 
increases in K g and X g . Thus, uncontrolled temperature 
variations can be responsible of part of the observed 
ageing of the physical properties of a granular material 
at rest. 

Our results also emphasize how delicate is the thermal 
cycling method compared to the classic tapping method 
1181 . even applied at low tapping amplitude lll9[ I2011 . 
In this sense, both methods we describe here certainly 
deserve further study and could be of great interest to 
probe the jamming transition of a granular assembly [9] 
and to unravel the local structure of grain displacement 
in the vicinity of this transition |[2lll . 

The mechanisms at stake in the compaction itself are 



still not well understood and deserve to be further inver- 
tigated. We already mentionned that the dilation of the 
container is not necessary, which indicates that the dif- 
ference between the thermal expansion coefficient of the 
beads and that of the container is probably not the pri- 
mary cause of the compaction, contrary to what was pre- 
viously proposed [8]. We are currently investigating the 
compaction of a granular column heated along its axis : 
A study of the response of the system as a function of the 
cycling period makes possible to analyze the role played 
by the penetration length of the temperature field (by the 
temperature gradient rather than by the temperature vari- 
ations) in the compaction process. 
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